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Abstract 

TeV-scale supersymmetry was originally introduced to solve the hierarchy problem 
and therefore fix the electrowcak (EW) scale in the presence of quantum corrections. 
Surprisingly, the numerical methods that evaluate the likelihood (or = — 21nL) 
to fit the experimental data and test the SUSY models do not account for fixing the 
EW scale. When this constraint is implemented, the likelihood (or x^) receives a 
significant correction (Sx^) that worsens the data fits of the models. We estimate this 
correction for the models: constrained MSSM (CMSSM), models with non- universal 
gaugino masses (NUGM) or higgs soft masses (NUHMl, NUHM2), the NMSSM and 
the general NMSSM (GNMSSM). Except the GNMSSM model and for a higgs mass 
rrih w 126 GeV, one finds that in these models 5x'^/ndf> 1.5, which violates the usual 
condition of a good fit already before fitting the observables, other than the EW scale 
itself {ndf=mxniheT of degrees of freedom). Given its significant (negative) implications 
for SUSY models, it is suggested that future data fits properly account for this effect, 
if one remains true to the original goal of SUSY. 
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1 A correction to the likelihood from fixing the EW scale. 



The main motivation for introducing low-energy (TeV-scale) supersymmetry (SUSY) was 
to solve the hierarchy problem and therefore fix the electroweak (EW) scale which remains 
stable under the addition of quantum corrections up to the Planck scale. This is achieved 
without a dramatic tuning of the parameters specific to the Standard Model (SM) [1]. 
SUSY is now under intense scrutiny at the LHC and the recent Higgs-like particle discovery 
of mass nifi near 126 GeV [2] brings valuable information on the physics beyond the SM. 
In this work we consider popular SUSY models and examine the impact on their precision 
data fits of the recently measured combined with the requirement of fixing the EW 
scale (vev v or mz) that SUSY was supposed to enforce, when it was introduced. 

Current likelihood methods to fit the experimental data and test SUSY models at the 
LHC involve a remarkable amount of technical expertise and work, backed by impressive 
computing power. Surprisingly, there is one important problem that these numerical meth- 
ods overlook when fitting the EW data. More exactly, the current likelihood to fit the data 
or its = — 21nL do not account for the "cost" of fixing the EW scale to its known 
value, that motivated the idea of supersymmetry. Using our previous [3], we show that this 
leads to an underestimate of the overall value of x^/^df in all popular models used at the 
LHC. Here n^f is the number of degrees of freedom of the model, defined as the number 
of observables fitted {no) minus that of the parameters of the model (up). The models 
analyzed include MSSM-like with different boundary conditions for gaugino and higgs soft 
masses, NMSSM and a generalized version of it, the so-called GNMSSM. 

For a realistic model, a good fit is expected to give a r atic0 

X^/ndf ~ 1, where Udf = no - n^. (1) 

In this work we show that in all SUSY models, the likelihood (or its x^) receives a 
correction due to the condition of fixing the EW scale to the measured value (m^), that 
worsens the current data fits. We show how to compute this correction (hereafter denoted 
(5x^) in the presence of the EW minimum constraints of the Higgs potential. The value 
of depends strongly on the value of the higgs mass m/j. Using the recent LHC result 
mil ~ 126 GeV, we find in the models other than the GNMSSM, a correction /ndj > 1.5 
without including the usual cost due to fixing observables other than the EW scale 
itself. Therefore, under the assumption of a simultaneous minimization of both 6x^ and the 
"usual" x^ (assumed to respect x^ jn^j ~ 1) these models have a total (x^ + '^X^)/'^d/ > 2.5, 
hardly compatible with the data (in the GNMSSM (x^ + ^X^)/™^/ ~ 2). 

Given its implications for these models and for SUSY in general, it is thus suggested 
that this correction be included in future data fits. In the following we substantiate these 
claims and analyze the consequences for the viability of popular SUSY models. 

1 For a set of observables Oi, one defines = Ei {Of - Of^fjoi. 
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2 The calculation of the correction to x^- 



Let us briefly review the numerical calculation of the likelihood in a SUSY model (see 
for example [H [5])- One usually chooses a set of observables Oi well measured such as: 
the W-boson mass, the effective leptonic weak mixing angle ^g^j, the total Z-boson decay 
width, the anomalous magnetic moment of the muon, the mass of the higgs {rrih), the dark 
matter relic density, the branching ratios from B-physics, Bs—Bg mass difference and also 
additional bounds (not shown) which should be counted in ridf^ too: 

mw, sin^ e'-^J'j, Tz, (^a^, nih, ftoMh^, 
BR{B -> Xs-y), BR(Bs fJ^^n'), BR{Bu tv), AMb,, etc. (2) 

Note that fixing the EW scale to its accurately measured value (m^) is not on this list, even 
though this motivated SUSY in the first place. For each observable Oi the corresponding 
probability is often taken a Gaussian P{Oi\^,y) where by 7 we denote the set of SUSY 
parameters that define the model, while by y we denote nuisance variables such as Yukawa 
(of top, bottom, etc) and other similar couplings. One then assumes that the observables 
are independent and multiplies their probability distributions to obtain a total distribution, 
which regarded as a function of 7, y (with "data" fixed), defines the likelihood L: 

L = Y[P{Oj\-f,y); 7 = {mo,mi/2,/xo,mo,^o,-Bo, • • • }; y = {yt,yb,- ■ ■} (3) 
j 

in a standard notation for the SUSY parameters, that are components of the set 7. To work 
with dimensionless parameters, all 7 should be "normalized" to some scale (e.g. the EW 
scale vev vq = 246 GeV). In the "frequentist" approach one maximizes L or equivalently 
minimizes the value of that, under a common assumption of Gaussian distributions, is 
defined as 

x' = -21nL. (4) 

One then seeks a good fit, such that x^/'^df ~ 1 at the minimum, by tuning 7, y to fit Oi. 
So, under the current view, one is tuning 7, y to fit Oi and this brings a contribution to 
the /iT-df "cost" in ([!]), but tuning the same 7, 1/ to fix the EW scale does not bring a 
contribution to /udf, it only brings a "fine" tuning! Such distinction made between the 
"costs" of these tunings is inconsisteno 

In other approaches like the "Bayesian" method one further combines L with priors 
(probabilities for 7, y) to obtain the posterior probability density [H El [71 [8]. In this method 
one searches for the point with the largest probability in parameter space, given the data. 

While the observables Oi are independent, the SUSY set of parameters 7 used to fit 
them, are not. They are usually constrained (correlated) by the EW minimum conditions. 

^ Also there is no technical difference or reason to treat separately these two "tunings" . 
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Indeed, in MSSM-like models, there are two minimum conditions of the scalar potential; 
one of them is determining the EW scale v as a, function of the parameters 7 of the model 
(it is not fixing v to any numerical value). In practice however, when doing data fits, one 
usually replaces v by hand, by the measured mass of Z boson (m^) and solves this minimum 
condition for a SUSY parameter instead (usually /Uq)- Doing so misses the impact on L of 
the distribution fixing the observable mz oc v to its measured m^, in the presence of the 
EW minimum constraints. The difference is important. The second minimum condition is 
fixing one parameter (say tan /3) as a function of the remaining 7. 

Let us then fix the electroweak scale in the presence of the two minimum conditions 
of the scalar potential for MSSM-like models and evaluate the total likelihood. To fix the 
notation, the potential is 

V = ml\Hi\'^ +ml\H2\'^ -(m^^i.i/s + h.c.) + {Xi/2)\Hi\^+ (A2/2)|/?2|^+ Agji^i 

+ Xi\Hi.H2\^ + [{X5/2){H^.H2f + Xg\Hi\\Hi.H2) + X^\H2WHi.H2) + h.c] (5) 



then denote 



A = Ai/2 cos^/3 + A2/2 sin^/3 + (A3 + A4 + A5)/4 sin2 2/3 + (Ae cos2/3 + A7 sin^/?); 
rn? = cos^ /3+ sin^ /? — m| sin2/3, (6) 

A is the effective quartic higgs coupling and m? is a combination of the higgs soft term^. 

To compute the overall likelihood, hereafter denoted Ly^{'y,y), that also accounts for 
the effect of fixing the EW scale, we use the results of |3] (see also [6]) that we adapt here 
to a form amenable to the calculation of the profile likelihooclfl. To find L^i'j, y) it is useful 
to write the two EW minimum conditions as Dirac delta of /i^2 defined below 

fi{T,y,v,P) ^ v-{-ni^/Xf'\ 

f2{T,y,v,/3) = tan/3 - tan/3o(7,y,v), (7) 

/3o is the root of the second minimum condition; /3o,/i,2 depend on the arguments shown; 
A, also depend on 7, y, /3. Taking account of constraints ([7]) the total ("constrained") 
likelihood L^(7,y) that also accounts for fixing the EW scale is 



Lw{l,y) = j dvd{iwiP) 5[fi{r,y,v,P)\ 6[f2{'y;y,v, f3)] 5{mz - m%) L{-f;y,v, /3) 
= vq L{r,y,vo,/3Q{'y,y)) S[fi{r,y,vo,/3o{j,y))] (8) 



^There is an extra condition in the case of NMSSM and GNMSSM models. 

*One can choose another parameter (instead of tan /3) like Bo, but the effect is just a change of variables. 

^ The EW min conditions are = -m'^/X, 2\{dm^ /dp)-m^{d\/dl3) = Q. See also [HE]. 

® The profile likelihood is found by maximizing the likelihood wrt nuisance parameters y, for fixed 7's. 
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Here vq = 246 GeV, m\ = gvQ/2 w 91.2 GeV, mz = gv/2, g^ = gl+g2 with 51^2 couplings 
for U(l), SU(2). To simplify notation, we did not display the numerical argument vq of /3o, 
i.e. we denoted /3o(7; v) = Pq{i-, U, vq). in rhs compensates the dimension of 6{mz—'m^). 

Li^j, y, V, /3) in the rhs of ([SD denotes the usual likelihood associated with fitting the 
observables other than the EW scale (mz), see eq.Q, before v, tan/3 are fixed by the 
EW min conditions. Integrating over v, tan/3 in the presence of the delta functions 
5(/2) is just a formal way to solve these EW minimum constraints and thus to eliminate 
these parameters in terms of the rest. We also introduced a 6{mz — m^) distribution, not 
imposed in data fits, which fixes the EW scale by enforcing the measured mass of Z boson 
(m^) and thus the replacement v — )• vq. Since is very well measured, using 5{mz — m^z) 
is indeed justified (however, the result can be generalized to a Gaussiar0). 

A comment about normalization in ([8|): one has to make sure that the two constraints 
(/i,2) are indeed normalized to unity wrt the variables over which we integrate them, in 
this case tan /3 and v, and above this condition is indeed respectecjl. Therefore we stress 
that there is no freedom for any additional factors to be present in eq.®. 

Eq.dSI) tells us that the original likelihood in the rhs, evaluated on the ground state 
and used in data fits is now multiplied by a Dirac delta evaluated at vq, 6{fi{...,vo, ...)). 
Further, one recalls that Dirac delta of a function fi{zi), i = 1,2, ...n, is related to that of 
its arguments according to a relation (formally derived by a Taylor expansion): 

6{Mz,)) = S[h{zf) + (V/i)o.(z - ^) + ■■■]= 6[n, {z, - z% (9) 

|V/i|o 

A summation over repeated index j is understood and the subscript "o" of |V/i|o means 
this quantity is evaluated at the point Zi = z^ where /i =0; rij are components of the normal 
n to the surface fi{z^) = 0, so n = (V/i/|V/i|)o. Expression ^ for 5 function together 
with eq.([8]) tell us that it is not enough for the parameters of the model to respect the 
constraint /i = which is what one actually checks in the numerical fits, and that there 
is instead an additional factor generated, represented by the gradient of the constrainllfl. 
The "constrained" likelihood of ([8]) becomes, after using ([9]): 

5rnj(lnzj — InzJ')] , , \\ 

Lwh,y) = ^ . 0^ L[r,y, vo,(3o{j,y)), Zi = {ij,yk} (10) 

Here 7^ , yk denote components of the sets 7 and y defined in ([3]) and we used ([9]) for In Zi 
instead of Zi as variables, to ensure dimensionless arguments for 5 function. denotes 
the absolute value of the gradient of /i evaluated at a point z^ = {7?, y^} that is a solution 



Additional distributions can also be considered for nuisance variables, like top, bottom masses, etc, 
and assumed to factorize out of L in the rhs, with corresponding integrals over the Yukawa couplings, see 
13]; however, these are not considered if one is interested to calculate the profile likelihood. 

*5(/i,2) are just distributions for i;,tan/3 (albeit special ones), so they also must be normalized to unity. 

® The gradient measures how the constraint "accommodates" the relative variations of the parameters. 
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to the EW min condition fi{'y ,y ,vq, j3o{'y ,y )) = 0, and has the value: 



^ \d\Q.Zi)o ■^V31n7o7o ^ \dlnyu Jo^ .. 

The subscript "o" stands for evaluation on the ground state (7^ = 7^^, yu = 1^1) ■ that 
emerged above has some resemblance to what is called the fine tuning measur wrt all 
parameters, both 7 and y. The arguments of Aq{'y^,y^) denote the parameters wrt which 
is computed and includes SUSY parameters and nuisance variables (Yukawa, etc). 

Eq. ()10p is the result in terms of distributions and can be written in an alternative form. 
The 6 in the rhs of (jlOp tells us that the Ihs is non-zero when Zi = for all i (i.e. 7^ =7^, 
yk = yk)- Another way to present this eq is to do a formal integration of (|lUp in the general 
direction Inz = rijlnzj, (sum over j understood, j running over the sets 7, y); this allows 
all independent parameters to vary simultaneousl'J^. After this, eq. (fTO]l becomes: 



LUi,y) = • (12) 



L in the rhs is exactly the usual, "old" likelihood computed in the data fits and evalu- 
ated at the EW minimum (refiected by its arguments) but without fixing the EW scale, 
while the Ihs does account for this effect. Note that due to the minimum condition 
fii'j^ ,y^ ,vo, /3o{'y^ ,y^)) = 0, one element of the set 7*^, say 7J?, becomes a function of 
the remaining, independent 7?, i ^ k. Usually 7^^ is taken to be /Uq of ([3]). 

The result in ec| . fjl2p shows that for a good fit in the model (i.e. maximal -L^^), one 
has to maximize not the usual likelihood in the rhs, but actually its ratio to A^. Let us 
introduce the notation Xw = — 21nL^ and also = — 21nL, then eq. (fT2]) becomes 

X^ = x' + 21nA,. (13) 

Therefore, after fixing the EW scale the usual receives a positive correction that 
depends on Aq and that is not included in the precision data fits. One can go further 
and compute from (|12p the profile likelihood by maximizing the likelihood (minimize Xw) 
wrt nuisance variables y^, for a fixed set of 7*^. This profile likelihood then takes the form 
L^(7°, y°(7°)) and can be used in detailed numerical analysis. 

Let us remark that in eq. (|10p we used In Zi as arguments under the Dirac delta function, 
which implicitly assumes that these are more fundamental parameters than Zi themselves, 
(here Zi = {"fj^yk})- In principle this is a choice, motivated here by the fact that it ensured 
dimensionless arguments for the delta function in pup. (unlike y^, Jj are dimensionful 

^"This was introduced in [lU] , with Amax = max^ |91n{j/91n7i|; see [TT] for an interpretation of 1/Amax- 
'^^Such integral is just a formal way of saying we solve /i = in favour of one particular 7I? (usually /lo). 
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parameters). Going from these parameters to their log's is a one-to-one change that does 
not affect the minimal value of totalil x^- If one insists in working with Zi as fundamental 
parameters, one simply changes In Zi — > Zi in eqs. ljlOp . (jlip . after "normalizing" 7, to some 
scale (e.g. vq), to ensure dimensionless arguments for 5 of (|10p . The result is that Aq is 
then computed wrt 7^ and yk instead of their logarithms, so in eqs. (fT2|) . (fT3]) one replaces 



Compared to its counterpart in Aq of (jlip . the second sum above is actually larger in this 
case since y/c < 1. In the following, for numerical estimates we shall work with Aq. A 
detailed investigation of the correction to is beyond the purpose of this paper and in 
the following we restrict the study to a numerical estimate. The main point is that 
receives a correction that needs to be taken into account. 

The above discussion can be extended to the Bayesian approach which is just a global 
version (in parameter space) of ()12p . In this case one assigns, in addition, some initial 
probabilities (priors) to the parameters of the model (7, y) then integrates over them the 
likelihood L.^^ multiplied by the priors. This gives the global probability of the model or 
"evidence", p(data), that must be maximized. The result is (see eq.(12) in [B] and Section 
2.2 in [3]): 

p(data) = / — - L(7,?/,z;o,^) X priors(7,y) (15) 

where the integral is over a surface in the parameter space 7, y defined by /i = /2 = 0. In 
this case 1/ Aq is itself an "emergent", naturalness prior, independent of and in addition to 
the original priors of the model. If one makes a change of a subset of variables to a more 
suitable one for calculations, an additional Jacobian factor emerges under the integral. We 
do not discuss further the Bayesian approach, for more details see [H El [TJ |8l [12l [131 El- 

3 A numerical estimate of the correction to SUSY. 

In this section we estimate in some models the correction dx^ obtained from eq. (|13p 

6xVndf = {2/ Hdf) In Aq (16) 

By demanding that the correction bx^ jn^] be small enough not to affect the current data 
fits results for x^ /ndj, one has a model-independent upper bound: 



^^In a Bayesian language, this would correspond to chosing log priors instead of flat ones for the param- 
eters. Ultimately this may reflect a problem of measure that is beyond the purpose of this work. 
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Ag < exp {ndf/2) 



(17) 



When this bound is reached, then 'X^/ridf = 1 + X^/?^4f) aiid assuming values of "usual" 
/n^f ~ 1 i.e. a good fit without fixing the EW scale, gives xt/fT'df ^ 2. 
Using p6|) . one can find the equivalent numerical value of that corresponds to one 
observable having a given deviation from the central value. This can be read below: 



3a 




A,p 


a 100. 


3.5 a 




A,p 


a 1000. 


5a 




A,p 


s 1 000 000 



(18) 

This gives, in a model independent way, a different perspective and a probabilistic inter- 
pretation to the numerical values of Ag (related to fine-tuning/naturalness) that avoid 
subjective criteria about this topic. As a result, any model with A^ > 100 (Ag > 200) 
is more than a 3a (3.25cj) away from fixing the EW scale. This is the case of the SUSY 
models that we discuss next. This effect is not accounted for by current fits. 

For our numerical estimates of 5%^ we restrict the analysis to using minimal values of 
Aq in SUSY models. Further, we only evaluate Aq wrt 7 parameters; notice that 

A,(7,y) > A,(7). (19) 

where the arguments display the variables wrt which Aq is actually computed. That is, 
we shall ignore the contribution to Aq due to variations wrt Yukawa couplings and other 
nuisance parameters. This underestimates the correctioJ^ Sx^- 

We consider the most popular SUSY models used for searches at the LHC, listed below: 

• the constrained MSSM model (CMSSM); this is the basic scenario, of parameterJ^ 
7 = {uiQ, m,i/2, fJ-o, Aq, i?o}, in a standard notation. Then Aq is that of ([TT]) with summation 
over these parameters only: 

For recent, two-loop numerical estimates of Aq we use the result in 

• the NUHMl model: this is a CMSSM-like model in which one relaxes the Higgs soft 
masses in the ultraviolet (uv), to allow values different from mg: = m^^ ^ ttt-q, with 
parameters 7 = {mo, mi/2, fJ^o, Aq, Bq, m™}. Aq is as in ([20]) with summation over this set. 
For two-loop estimates of Aq as a function of the two-loop leading-log ruh we use [6] . 



Actually, if one integrates (fixes) the masses of the top, bottom, etc, to their numerical values, it is 
Aq{'y) that is finally obtained in ea. (|13|) Indeed, integrating in (O over additional d(lnj/t) and d(lnyi,) 
and with distributions S{mt — to?) and 5{mi, — to?) to "fix" the accurately measured observables, one obtains 
a result similar to H10|) but no contribution is present in egs. HllI) to H13[) from Yukawa couplings [3]. 

^^Parameters 7 are those following from ()10p . and this is why /io is quoted instead of usual sgn{^o). 
Similar for the other models. Also Bo is quoted instead of tan/3, the difference is a change of variables. 
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• the NUHM2 model: this is a CMSSM-hke model with non-universal Higgs soft masses, 
m]^^ T^m^^ T^mo, with independent parameters 7 = {mQ,mi/2, fJ-o, ^o,Bq, m™, m^^}. Then 
Aq is that of (j20p with summation over this set. For two-loop numerical estimates of Ag 
we use the results in ref. [B]. 

• the NUGM model: this is a CMSSM-like model with non-universal gaugino masses mx- , 
i = 1,2,3, with 7 = {niQ, fiQ, Aq, BQ,mx-^, m\2,m\^}. Ag is that of ([20]) with summation 
over this set. For two-loop numerical estimates of Ag we use the result in [6]. 

• the NUGMd model: this is a particular case of the NUGM model with a specific relation 
among the gaugino masses m\-, i = 1,2,3, of the type = r]imi/2, where take 
only discrete, fixed values. Such relations can exist due to some GUT symmetries, like 
SU(5), SO(IO), etc [15]. The particular relation we consider is a benchmark point with 

= (l/3)mi/2) "^Ai = (— 5/3)mi/2; rnx^ = rni/2, corresponding to a particular GUT 
(SU(5)) model, see Table 2 in [15]. is that of ([20]) with 7 = {mo, mi/2) ^O) -Bo, /^o}- For 
numerical estimates of we use ref. [6]. 

• the next to minimal MSSM model (NMSSM): the model has an additional gauge singlet. 
Its parameters are 7 = {mo, fiQ, Aq, BQ,mi/2i^s}i with ms the singlet soft mass. For an 
estimate of A^ we use the results in [16l |T7] . Notice that these papers evaluate instead 
Amax which is the largest fine tuning wrt to any of these parameters (instead of their sum 
in "quadrature" as in Ag). As a result, Amax is usually slightly smaller, by a factor between 
1 and 2 as noticed for the other models listed above [6] and thus dx^ is underestimated. 

• the general NMSSM (GNMSSM) model: this is an extension of the NMSSM in the sense 
that it contains a bilinear term in the superpotential for the gauge singlet, MS'^, see for 
example |16j . So the singlet is massive at the supersymmetric leveO and we have an 
additional parameter M to the set we have for the NMSSM. Again, A^ax is used here [16] 
instead of Ag, so 6x^ is again underestimated (as for the NMSSM). 

Our estimates for Sx^/ndf for different higgs mass values are shown in Tables [1] and [2j 
The results present the minimal value of Ag evaluated in the above models as a function of 
the higgs mas^il, after a scan over the entire parameter space (all 7, y and also tan (3 of the 
corresponding model, as described i [6]). We show the values of Ag for central values 
of nih close to the pre-LHC lower bound ~ 115 GeV in Table [Tj and for 123 to 127 GeV in 
Table [2l These values allow us to account for a 2-3 GeV error of the theoretical calculation 
at 2-loop leading log level |19j . For rrih ~ 115 GeV, one could still have Sx^/n^f < 1 for 
some models, that could have allowed a corrected x^/ndf ~ 1. Further, minimal Ag grows 
approximately exponentially wrt to nih, due to quantum corrections. Indeed, Ag ~ rn^usy 
and since the loop correction is roughly 5mh ~ Inmsusy one finds Ag k, exp(m/j/GeV). As 
a result, a strong variation of Ag wrt is found, and Sy^ increases by w 1 for a 1 GeV 
increase of mh, see Table [2j 

'^^Its mass can be of few (5-8) TeV, so one can integrate it out and work near the decoupling limit |18| . 
^® See Figures 1 to 8 in [6] showing as a function of vrih, after scanning the entire parameter space. 
^^The scan included a range of [—7, 7] TeV for Aq, mo and mi/2 up to 5 TeV, and 2 < tan/3 < 62 and 
also allowed a 2a deviation for fitted observables and for a 3a for Q.oMh^ , see for details [Bj. 
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Model 




Approx 


y 


5X^(115) 




CMSSM 


5 


2-loop 


15 


5.42 


9 


NUHMl 


6 


2-loop 


100 


9.21 


8 


NUHM2 


7 


2-loop 


85 


8.89 


7 


NUGM 


7 


2-loop 


15 


5.42 


7 


NUGMd 


5 


2-loop 


12 


4.97 


9 


NMSSM 


6 


1-loop 


12 


4.97 


8 


GNMSSM 


7 


1-loop 


12 


4.97 


7 



Table 1: The correction 6x^ = 21nAq(7), the number of parameters Up and degrees of freedom 
Udf in SUSY models, for irifi w 115 GeV corresponding to the pre-LHC bound of m/j. Notice that 
in this case Sx^/ndf < 1. ridf may vary, depending on the exact number of observables fitted. The 
numerical values of A, are from [B] (first 5 models) and [H] for NMSSM, GNMSSM. 



Model 


A. 


(5x^(123); 




5X^(125); 


A, 


5X^(126); 


A, 




CMSSM 


380 


11.88 


1100 


14.01 


1800 


14.99 


3100 


16.08 


NUHMl 


500 


12.43 


1000 


13.82 


1500 


14.63 


2100 


15.29 


NUHM2 


470 


12.31 


1000 


13.82 


1300 


14.34 


2000 


15.20 


NUGM 


230 


10.88 


700 


13.10 


1000 


13.82 


1300 


14.34 


NUGMd 


200 


10.59 


530 


12.55 


850 


13.49 


1300 


14.34 


NMSSM 


>100 


9.21 


>200 


10.59 


>200 


10.59 


>200 


10.59 


GNMSSM 


22 


6.18 


25 


6.43 


27 


6.59 


31 


6.87 



Table 2: As for Table 1, with and corresponding (5%^ = 21nAq(7) for rrih equal to 123, 125, 
126 and 127 GeV (shown within brackets). This can also show the impact of the 2-3 GeV error 
in the theoretical calculation of nih [HI HO] • grows « exponentially with rrih [3 [S] ; a 1 GcV 
increase of m/i induces about 1 unit increase of Sx^- In all cases except GNMSSM, fixing the EW 
scale brings a correction by^ jridf > 1. Note bx^ judf can be larger if one also includes the impact of 
Yukawa couplings on Ag. Same loop approximation, number of parameters and degrees of freedom 
apply as in Tabled! The values of A, are from [B] and from [H] for NMSSM, GNMSSM. 

In all models, for the currently measured mh ~ 126 GeV, 5x^ /udf alone is larger than 
unity (or close to 1 for GNMSSM), without considering the original contribution due to 
the "usual" 'X^ jn^^ coming from fitting observables other than the EW scale. The above 
correction is too large by the usual criteria that total x^/udf ~ 1. Further, assume that in 
the above models one finds a point in the parameter space for which x^/'^df ~ 1 and then 
add to it the effect of Sy^. One then finds that X^/n^f is close to or larger than 2, see 
Table O In this Table, the values of reduced x? increase (decrease) by ~ 0.1 for an increase 
(decrease) of nih by 1 GeV, respectively, except in the GNMSSM where this is even smaller 
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mh = 126 GeV: 


CMSSM 


NUHMl 


NUHM2 


NUGM 


NUGMd 


GNMSSM 


Xl/ndf ■■ 


2.66; 


2.83; 


3.05; 


2.97; 


2.49; 


1.94 



Table 3: An estimate for total Xw/^df in various SUSY models for m/i«126 GeV and x^/n^/wl. 



(0.04 for 1 Gev). The log dependence (5x^ = 21nAq means that uncertainties in evaluating 
Aq are reduced and 6x^ values listed in Table [2] and used in Table [3] are comparable, even 
though corresponding Ag's are very different. 

Our estimates for Xw/''^df shown in Tableware hardly acceptable for a good fii^. Inter- 
estingly, increasing the number of parameters of a model (decrease ridf) could decrease 
and (5x^, but this reduction may not be enough to reduce Xw/''^df, since ndj is now smaller. 
This is seen by comparing the reduced i^i NUGM and CMSSM. This is because 6x^ 
depends only mildly (log-like) on A^, so only a significant reduction of Aq can compensate 
the effect of simultaneously reducing n^j. For such case compare CMSSM with GNMSSM. 

The values of Xw/''^df could be higher than our estimates above, since they ignore that: 
a) we used only minimal values of A^ over the whole parameter space, b) Yukawa effects 
on dx^ were ignored and these can be significani^. There is also an argument in favour of 
a smaller Xw/^df^ that current theoretical calculations of nih may have a 2-3 GeV error. 
Assuming this, for rrih ~ 123 GeV (instead of 126 GeV), in GNMSSM one obtains a small 
change: Xw/''^df = 1-8 while for the other models this ratio is >2.3. Another reduction 
may emerge in numerical analysis if using eq. (|14p instead of (jlip . but the impact of its 
larger Yukawa contributions (that is usually also larger than that due to 7's) makes this 
possibility less likely. 

Regarding the NMSSM model. Table [2] only provided a lower bound on Aq. However 
we can do a more accurate estimate, using recent data fitj^ that evaluated both x^ ai^d 
Aq. To this purpose, we use the minimal value for x'^ = 6.4 in [T7j (last two columns of 
their Table 3) together with its corresponding^ Ag = 455. This gives a x^/udf = 18.64/8, 
which is similar to the other models discussed abovJ^. Further, using instead our Table [2] 
where 5x^ > 10.59 we find Xw/^df > 2.32 (for x'^ /^df ^ 1); which is in agreement with the 
aforementioned value derived from accurate data fits. 

To conclude, the requirement of fixing the EW scale has a strong effect on the value 
of Xw I ''^df ) with negative impact on the data fits and on the phenomenological viability of 
these models. While our numerical results are just an estimate of the correction Sx^, the 
effect is nevertheless present and demands a careful investigation by the precision data fits. 

Assuming a distribution, the p- value in these models would be < 1% (and 5% for GNMSSM). 
^® In CMSSM, Ag wrt top Yukawa alone is larger than that wrt to any SUSY parameter, see fig. 2 in [5]. 
^°For other recent data fits see [2ll [22l [23l [24] . 

Aq is even higher as it is not computed according to (|20|) but reports max values wrt each parameter. 

Other values quoted in [T7] bring an even larger value for this ratio. 
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4 Some implications for model building. 



Let us discuss some implications of the above result for model building. 

a) . A natural question is how to reduce /ridf- Here are three ways to attempt this: i) 
additional supersymmetric terms in the model, which unlike SUSY breaking ones, are less 
restricted by experimental bounds; ii) additional gauge symmetry, iii) additional massive 
states coupled to the higgs sector. All these directions have in common a possible increase 
of the effective quartic higgs coupling (A) so one can more easily satisfy an EW minimum 
conditiorl^ v"^ = —m?/\ for v ~ O(lOOGeV), m ~ O(lTeV), that demands a larger A. As 
a result one can obtain a smaller Sy^ oc In Ag. As mentioned, the increased complexity of 
the model (more parameters) is to be avoided, since then Udj can decrease and 5x^ /udf 
may not change much (or even increase). The GNMSSM model is an example of i), where a 
supersymmetric mass term for the additional singlet essentially enabled a smaller by^ Jn^f 
than in other MSSM-like modelS Similar but milder effects exist in the NMSSIvO- An 
increase of A could also be generated by using idea ii) by adding more gauge symmetry 
(e.g. [26] )• Regarding option iii), one can consider additional massive states that couple 
to the higgs sector, and that in the low energy generate corrections to the higgs potential 
and effective A and m/, [27] with similar effects. These ideas may indicate the direction for 
SUSY model building. 

b) . The above negative implications for some SUSY models remind us about the real 
possibility that no sign of TeV-scale supersymmetry may be found at the LHC. If so, this 
can suggest its scale is significantly larger than few TeV. Alternatively, one could attempt 
to forbid the existence of asymptotic supersymmetric states while trying to preserve some 
of the nice advantages of SUSY. In such scenario, superpartners would be present only 
as internal lines in loop diagrams, they would not be real, asymptotic final states. One 
could describe this situation by using some variant of non-linear supersymmetry that can 
be described in a superfield formalism endowed with constraints (see examples in |28|). 
The hope would be to preserve SUSY results like fixing the EW scale, gauge couplings 
unification, radiative EW symmetry breaking, for which the superpartners in the loops 
play a crucial rolJ^. However, it is difficult to realize this idea in practice. One also 
recalls the supersymmetric quantum mechanics case where SUSY is used only as a tool for 
performing complex calculations [30], which could suggest ideas for the field theory case. 

A tension in the relation v'^ = —rn^/X translates into a larger Sx^ ~ In Aq [S]. Recall that in MSSM A 
is very small and fixed by gauge interactions (at tree level) and this is one source for the above problems. 

^''For a large value of the supersymmetric mass term of the singlet (M of few TeV, 5-8), the correction 
to the higgs effective quartic coupling A is ~ (2/i/M) sin 2/3 already at tree level [TH], with impact on A, 
and Sx^- [IS]. For a recent study of the GNMSSM and its LHC signatures see |25| . 

The correction to A and nih is in this case restricted by perturbativity in the singlet coupling A (of 
\SH1.H2) and also its proportionality to Asin^ 2/3 instead, while in the GNMSSM is oc sin 2/3. 

related idea exists [29], based on a similarity to gauge fixing in gauge theories and the subsequent 
emergence of the ghost degrees of freedom of different statistics. This would have an analogue in the above 
SUSY scenario in "fixing the gauge" in the Grassmann space. Similar to the emergence of ghosts in gauge 
theories as non-asymptotic states one could attempt to obtain non- asymptotic superpartners states. 
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c). The remaining possibility is to abandon (low-energy) SUSY and eventually consider a 
different symmetry instead. One can use model building based on SM extended with the 
(classical) scale symmetry, thus forbidding a tree level higgs mass. This symmetry is broken 
at the loop level by anomalous dimensions, which would bring in only log-like dependence 
on the mass scales |3T]. Additional requirements (unitarity, etc) could be added. For model 
building along this direction see some examples in [32] and references therein. 

5 Conclusions 

The main motivation for TeV-scale supersymmetry was to solve the hierarchy problem 
and therefore fix the electroweak scale (vev v or mz) in the presence of the quantum 
corrections. Rather surprisingly, the numerical methods that evaluate the likelihood (or 
its = — 21nL) to fit the data in SUSY models do not account for the "cost" that is 
due to fixing the EW scale to its measured value {m^). When this condition is imposed 
and combined with the constraints of the EW minimum conditions of the higgs sector, one 
finds that receives a positive correction, 5x^ > 0. This correction should be included in 
the analysis of the total of the supersymmetric models. 

This result is important, since data fits report separately the result for minimal x^/f^d/ 
that is a measure of tuning the parameters to fit some observables (other than the EW 
scale), and the amount of another tuning (a 'fine" one)0 needed to fix the EW scale, 
whose x^ "cost" is ignored. This suggests an inconsistency. Our results show that these 
two aspects ("tunings") are actually part of the same problem rather than separate issues 
and thus the total x^ + ^X^ must be examined by the data fits of the models. 

For the recently measured value of the higgs mass (~ 126 GeV), the correction dx^ 
was estimated and was shown to be significant in most popular SUSY models: constrained 
MSSM (CMSSM), models with non-universal higgs soft masses (NUHMl, NUHM2) or 
with non- universal gaugino masses (NUGM), in the NMSSM and its generalized version 
GNMSSM. This correction has negative implications for the data fits of SUSY models. 
Our estimates show that for ruh ~ 126 GeV, this correction alone is ^x^ jn^^ > 1.5, which 
violates the traditional condition for a good fit already before fitting observables other 
than the EW scale. Adding this to a good fit according to the "usual" x^ /ndf ~ 1 would 
give (x^ + '^X^)/^d/ > 2.5, hardly acceptable. Further contributions to 5x^ also exist from 
Yukawa couplings, not discussed. The effect is however milder for the GNMSSM case. 

To conclude, the requirement of fixing the EW scale in SUSY models brings a cor- 
rection to the likelihood, with negative implications for the overall data fit and for the 
phenomenological viability of these models. While we only provided a numerical estimate 
of it, that needs further investigation, this correction is present and should be included in 
the total x^, if one remains true to the original goal of SUSY of fixing the EW scale. 

"measured" instead by various model dependent definitions. 
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